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Abstract 

We show that moderate deviations from the Maxwell-Boltzmann energy distribution 
can increase deuterium reaction rates enough to contribute to the heating of Jupiter. 
These deviations are compatible with the violation of extensivity expected from 
temperature and density conditions inside Jupiter. 



1 Introduction 

Jupiter emits more radiation then it recives from the sun: the origin of this 
excess heat is still uncertain and debated. Possible explanations are: release 
of gravitational potential energy due to the planet contraction and/or Helium 
sedimentation [1], decaying of radioactive isotopes in the core [2], or deuterium 
burning [3] . Each of these hypotheses has difficulties [4] ; in particular, standard 
calculations of deuterium burning reaction rates predict negligible contribution 
to the planet thermal balance, in spite of the substantial enhancement due to 
electron and ion screening effects [3,5]. 

In a strongly coupled plasma, anomalous diffusion and time correlation ef- 
fects originate non-Maxwellian two-body relative energy distribution that can 
be parameterized with and, in same cased, assume the same functional form 
that appears in the contest of Tsallis non-extensive thermodynamics [6,7]. As 
demonstrated for the solar core, small changes of the tail of the energy dis- 
tribution can strongly modify the fusion rates without affecting mechanical 
properties (hydrostatic equilibrium and the sound speed) that depend on the 
mean value of the distribution [8]. 

Since the internal conditions of Jupiter indicate the existence of a strongly 
coupled plasma, we investigate the effects of the consequent small deviations 
from the standard Maxwell-Boltzmann (MB) statistics on deuterium burning 
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rates and the possibility that deuterium burning could play, or have played in 
the past, a role in Jupiter thermal balance. 



2 Jupiter interior and standard deuterium burning 

Jupiter interior is a mixture of liquid metallic hydrogen and helium, with 
density and temperature within the ranges: p.j = 2 — 5 g cm -3 and kT =1 — 2 
eV; during planet formation the central temperature should have been of the 
order of 10 — 20 eV [4]. If we assume the reference values pj = 5 g cm -3 
and kT = 2 eV, the corresponding density of H, D, He and electrons are: 
n p = 2.4xl0 24 cm" 3 , n D = 3xl0 _5 n p , n He = 6.25xl0~ 2 n p , andn e = 1.125 n p . 

Being the electron Fermy energy, Ep ~ 50 eV , much greater than the thermal 
energy, the electron gas is fully degenerated. The electron and ion plasma 
parameters, Tj = (ZiZje^/^kT), where a = (3Z/(4vm)) 1/3 is the Wigner- 
Seitz radius, are: r e pa Tj pa 16. Therefore, the interior of Jupiter is a Strong 
Coupled Degenerated Plasma. 

Jupiter excess energy flux $ j = 5.4 x 10 3 erg cm -2 sec -1 implies an excess 
luminosity: Lj = 3.5 x 10 24 erg sec -1 . Each D(p,7) 3 He reaction produces 
an energy Q p p> = 5.493 MeV. Therefore, deuterium burning is relevant to 
Jupiter thermal balance only if the rate is greater than a threshold r t = 
(pjLj)/(MjQ pD ) pa lsec^ 1 cm~ 3 , where we have used as Jupiter mass and 
radius: Mj = 1.90 x 10 33 g and Rj = 7.14 x 10 9 cm. 

Two-body reaction rates in a thermal plasma can be expressed as (l+<5jj)rjj = 
Tbjjij ((tv) Mi where {<jv)m is the thermal average of the reaction cross section 
a times the relative velocity v, and rii is the number density of species i. 

Cross sections between charged particles below threshold are dominated by 
the penetration factor e~ b /^, where b = -n^plpZiZje 1 /% and p is the particle 
reduced mass. In fact the astrophysical factor S(E) = a(E)Eexp {b\J~E) has 
a mild dependence on E in absence of resonances. 

There exist no experimental determinations of S(E) for deuterium reactions 
at energies in the eV range. Theoretically motivated low-energy extrapolations 
can increase the deuterium [9] are: S(0) = 2.0 x 10~ 7 MeV b for D(p, 7) 3 He, 
S(0) = 5.0 x 10- 2 MeV b for D(d, n) 3 He and S(0) = 5.6 x 10~ 2 MeV b for 
D(d, p)T. 

If the thermal distribution is Maxwellian, (<jv)m ~ Jo°° Eo'(E)e~ E ^ kT dE, ther- 
monuclear rates are dominated by the product of the exponentials e~ E l kT (MB 
distribution) and e^ 6 /^ (penetration factor), product that has its maximum 
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at the Gamow Peak Energy, Eq. In the weak screening regime E s < Eq we 
can use a saddle point expansion around E and find: 



where E = (bkT/2) 2 ^ 3 is the Gamow Energy, J max = exp (— 3E /kT) is the 
integrand at E — E , and A = 4(E kT/3) 1 ^ 2 measures the width of the 
Gamow Peak. The two correction factors, f e — (1 — E s /E )exp (E s /kT) for 
electrons and /; = exp (H(0)/kT) for ions, take care of screening. 

In fact screening is very important in Jupiter interior. Electron screening low- 
ers the coulomb barrier of E s = e 2 /D s , where D s is the short-range screening 
distance, which is determinated through the incipient Rudberg method [5]. 
Many body correlations modify also ion wave functions in plasmas: it is pos- 
sible to define an equivalent mean field potential H(r) which is related to the 
plasma parameter T in the contest of Ion Sphere Model: H(0)/kT = 1.057r s 
where T s is the screened ion plasma factor T s = Ti exp (—a/D s ) and Tj is the 
ion plasma factor. 

Despite the large enhancement factors (f e ~ 10 11 , j\ ~ 10 2 ) deuterium burning 
remains negligible: r v r> ~ 1CT 19 and tdd ~ 1CT 26 sec -1 cm" 3 . Rates are much 
lower than the threshold rate r t for any relevant temperature. 



3 Non-Maxwellian deuterium burning 

Non extensive Thermodynamics, introduced by Tsallis [10] has been applied to 
many different fields. In particular, it has been observed [6,7] that small devia- 
tions from Maxwellian distribution can be described by weakly non-extensive 
Tsallis distributions, and that such deviations can have dramatic consequences 
for nuclear reaction rates in the solar core [8]. 

We can calculate the effects of non-extensivity on nuclear reaction rates by sub- 
stituting the Maxwell-Boltzmann distribution, exp (—E/kT), with the Tsallis 
one, [1 — (1 — q)E/kT] q ^ 1 ~ g \ inside the thermal average: 

7 b / £\«/(H) 
Mm ~ J e VWm q)—j S(E + E s ) dE . (2) 

o 

The limit q — > 1 recovers the Maxwellian case. Analogously to the Maxwellian 
calculation, we find a g-dependent Gamow Peak Energy E 0q , and, in the weak 
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screening condition E s < E 0q , we can compute again the reaction rates at first 
order in the saddle point expansion: 



n inj [2( E s \ S(E 0q ) 
l + ^VA Ej (A;T) 3/ 2 W^n 



(3) 



where now E 0q = E £ 2 (T, q)/q 2 ^ is the g-dependent Gamow Peak Energy, 



L g,max 



cxp 



^E 0q + E s 1-9 V 



is the integrand at E — E 0q , and 



A q = a far, q),^ 12 (i + gO^MTOW!^ 

V V ^ ' V (l-(l-q)(E 0q -E s )/kTf J 



-1/2 



measures the peak width. The other terms are: a = (1 — q)b 2 ^ 3 / (4g 2 A;T) 1 ^ 3 , 
6 = 1 + (1 - q)E s /kT, and 



3f (T, g) = - a + 2~ 1/3 (276 - 2a 3 - c) 3 + (276 - 2a 3 + c 



with c = ^276(276 -4a 3 ), while E and A are the same as in the MB case. 

Note that it is not possible to factorize an electron enhancement factor as for 
the MB distribution. 



4 Discussion and conclusions 



In Fig. 1 deuterium reaction rates are plotted for two different values of the q 
parameter. We can observe that moderate deviations (q — 1) ~ 0.1 from the 
Maxwellian distribution increase deuterium burning rates above the threshold 
rate r t as lsec _1 cm -3 : therefore, these processes should contribute to heat 
Jupiter at the present epoch (T 1 — 2 eV). If we read the graphs for T 10 — 
20 eV, which corresponds to temperatures of the planet during its formation, 
we realize that it was sufficient a smaller value, (q — 1) ~ 0.03, to effect the 
thermal balance in that period. 



Reaction rates of the order of that required to heat Jupiter do not cause 
decrease significantly the deuterium density inside the planet; in fact a burning 
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rate equal to ten times the threshold rate would consume a significant fraction 
of deuterium only after a time of the order of nc/(10r t ) ~ 10 11 years. 

These considerations demonstrate that deuterium burning is a possible expla- 
nation for the Jupiter excess heat. Precise determinations of the conditions 
inside Jupiter and additional microscopic calculations are necessary for a bet- 
ter determination of the range of values of q relevant to Jupiter interior. 
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q = 1 .03 





Fig. 1. Reaction rates for the reactions p+D (solid), D+D (dashed) and p+p 
(dot-dashed) as function of temperature for two values of the Tsallis parameter: 
q = 1.03 (upper frame) and q = 1.1 (lower frame). 
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